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THE SPECTRAL EVOLUTION OF THE FIRST GALAXIES. I. 
JWST DETECTION LIMITS AND COLOUR CRITERIA FOR POPULATION III GALAXIES 

Erik Zackrisson^'^*, Claes-Erik Rydberg^'^, Daniel Schaerer''"', Goran Ostlin^'^, Manan Tuli^'^ 

ABSTRACT 

The James Webb Space Telescope (JWST) is expected to revolutionize our understanding of the 
high-redshift Universe, and may be able to test the prediction that the first, chemically pristine 
(population HI) stars formed with very high characteristic masses. Since isolated population HI 
stars are likely to be beyond the reach of JWST, small population HI galaxies may offer the best 
prospects of directly probing the properties of metal-free stars. Here, we present Yggdrasil, a new 
spectral synthesis code geared towards the first galaxies. Using this model, we explore the JWST 
imaging detection limits for population HI galaxies and investigate to what extent such objects may 
be identified based on their JWST colours. We predict that JWST should be able to detect population 
HI galaxies with stellar population masses as low as ~ 10^ Mq at z « 10 in ultra deep exposures. 
Over limited redshift intervals, it may also be possible to use colour criteria to select population HI 
galaxy candidates for follow-up spectroscopy. The colours of young population HI galaxies dominated 
by direct star light can be used to probe the stellar initial mass function (IMF), but this requires 
almost complete leakage of ionizing photons into the intcrgalactic medium. The colours of objects 
dominated by nebular emission show no corresponding IMF sensitivity. We also note that a clean 
selection of population HI galaxies at z « 7 — 8 can be achieved by adding two JWST/MIRI filters to 
the JWST/NIRCam filter sets usually discussed in the context of JWST ultra deep fields. 
Subject headings: stars: Population HI - galaxies: high-redshift - cosmology: dark ages, reionization, 
first stars 



1. INTRODUCTION 

Due to the lack of efficient coolants in metal- 
free, primordial gas, the first stars to form in our 
Universe were likely very massive (> 10 Mp, and 
possibly even up to 1000 M(d\ e.g. Brom m et al.l 



19991: iNakamura &: Ume mura'12 0021: iTan fc McKeel l2004 
Greif fc Bromml l2006t lOhkubo etai.1 l2009f ). Some of 



these population HI stars probably exploded as super- 
novae (Hcger ct al. 2002) and initiated the chemical en- 
richment which allowed less massive and more metal- 
rich stars to form (the population I and II stars known 
from the local Universe, with characteristic masses below 
1 Mq). The highly energetic radiation emitted from pop- 
ulation HI stars may have played an impo rtant role in the 
cosm i c reio nization at rcdshifts z > 6 fe.g. Sokasia n et al.l 
12004 ITrenti fc Stiavelli 2009 ). and the remnants left be- 
hind by these objects may also h ave acted as seeds for the 
supermassive black holes (e.g. ITrenti fc Stiavellil I2007D 
known to exist already at z 6 (e.g. lWillott et al 

Current predictions for the initial mass function of pop- 
ulation HI stars are, however, still shaky. The most 
recent simulations suggest that the typical population 
HI masses may be lower than previously believed (e . ; 



IStacv et al.l[200l IClark et al.l[20Tlt IGreif et al.ll20TT£ 

and the chemical abundance patterns of extremely metal- 
poor stars in the Milky Way halo also point towards early 
supernovae with progenitor masses closer to 10 Mq 
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than ^ 100 AI^ (for rev iews, see iBeers fc Christliebl 
120051: IKarlsson et all 120 111 ). At the same time, it is 
not clear that the chemical signatures of the most mas- 
sive population HI supernovae would be detectable in 
current halo samples. The stars beari ng the mark of 
such s upernovae may be t oo scarce (e.g. ISalvadori et al.l 
[20071: ITrenti fc Sh{d|[20l 0). hiding at too high metal- 
licities (Karlsson et al.ll2008i) or i n the inner regions of 
the Galaxy fe.g. lTumlinsonll2010D . Direct observations 
of massive population III stars in the high-redshift Uni- 
verse would help settle the issue, but this is beyond the 
capabilities of current telescopes. 

The very first population HI stars are expected to 
form in isolation or in small numbers within ^'lO^- 
10^ M(7) dark matter halq s at redshifts z » 20-50 
(e.g. iTegmark et"an 119971: lYoshida et all I2003D . but 
the prospects of detecting suc h stars on an indi- 
vidual basis app e ar bleak (e.g. [Gardner et al~1 |2006| ; 
IGreif et all [20091 : iRvdberg. Zackrisson fc Sco tt 20l(f] 
at le as t before they go supern ovae (jWeinmann fc Lillvl 
[200l IWhalen fc Fr^ [20Tol ). However, popula- 
tion III stars may continue to form within the 
more massive halos (> 10^~* -^0) hosting the 
first galaxiej^ at z < 15 (IScannapieco et al.l [2(j03l: 
Schneider et al.l I2006al: i Tornatore. Ferrara fc Schneider 



20071: I Johnson. Gre if fc: Br omml 12008): [Johnson et al 



2009t iStiavelh fc Trentil I20i0l: lJohnson[ |2010'). and this 



could in principle allow their integrated signatures to 
be detected wit h the upcoming James We bb Space Tele- 
scop^ (Bromm et al.ll200"ll : [Johnson et al.ii2009, : Johnsonl 

^ Due to their low expected stellar population masses, these 
objects are sometimes also - and perhaps more appropriately - 
referred to as population III star clusters 

^ http://www.jwst.nasa.gov/ 
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doTp). 

Since the first generation of galaxies is predicted to 
form in high-density regions that have been pre-enriched 
by pop III stars in minihalos, these galaxies are not ex- 
pected to be metal-free, and ar e most lik ely dominated 
by chemically enriched stars (e.g. lGreif et al . 2010). True 
population III galaxies may, however, form at slightly 
later epochs, in low-density environment s which have re- 
mained chemically pristine (e.g.[S cannap ieco et al1l2003l : 
iTrenti et al.l[200a IStiaJelh fc Trentiil20ir 

Due to the exceedingly high temperatures of mas- 
sive population III stars {Tcs ^ 10^ K), such ob- 
jects may contribute substantially to the photoioniza- 
tion of the interstellar medium and alter the overall 
spectra of their host galaxies. A number of authors 
have pointed out that population Ill-dominated galax- 
ies could be identified based on the strengths of the 
Lya line, the HcII A1640, A4686 lines, or the Lyman 
'bump' (e.g. iT umlinson & Shull 2000; Tumlinso n et al 



2001 : Oh et al.li2001i.Malhotra fc Rhoads.2002^ .Schaerer 
200a 120031: llnoud 120101) " and several searches have since 
tried to apply these techniques (e.g. 'Dawson et al .ll2004l: 
Nagao et al. 2005;^ Dijkstra fc Wyith e 2007; Ouc hi et al l 
2008t [Nagao et al.l 120081: ICai et al.l 1201 iD . While a 
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exact nature remain unclear. So far, essentially all 
searches for population III galaxies have relied on spec- 
troscopy or narrowband photometry. If population III 
galaxy candidates could instead be singled out from large 
samples based just on their broadband colours (as re- 
cently attempted bv lBouwens et al.ll2010D . this would al- 
low for substantial gains in terms of observing time. 

Here, we present Yggdrasil, a new spectral synthesis 
model for the first generations of galaxies, and use it 
to explore the spectral signatures of pure population III 
galaxies in James Webb Space Telescope (JWST) photo- 
metric surveys. The technical details of the model are 
described in Sect. [2] In Sect. |3l we define three differ- 
ent classes of first galaxies, which differ by the relative 
amount that nebular emission is expected to contribute 
to their integrated spectra. The resulting James Webb 
Space Telescope (JWST) mass detection limits for both 
population III and more chemically evolved galaxies in 
photometric surveys are presented in Sect. ID Potential 
strategies for searching for population III galaxies (dom- 
inated either by nebular emission or direct star light) in 
deep JWST images are presented in Sect. [5] and Sect. [SI 
A number of potential caveats with the proposed ap- 
proaches are discussed in Sect. [7] Sect. [8| summarizes 
our findings. 

2. THE YGGDRASIL MODEL 

While the properties of galaxies at 2 > 10 remain 
an uncharted territory observationally, theory and sim- 
ulations provide a number of clues as to what one 
might expect from such objects. In the models of 
[Stiavelli fc 'Ren ti (2010), the very first galaxies are pre- 
dicted to form in high-density regions that have been 
pre-enriched by population III stars that exploded as su- 
pernovae in minihalos at even higher redshifts. True 
population III galaxies start to form somewhat later, 
in low-density environments which have remained chem- 



ically pristine. Since pockets of primordial gas may 
survive in galaxies that have already experienced some 
chemical enrichment, hybrid galaxies in which popula- 
tion III, population II (representative of the metallicity 
in the Milky Way halo; Z < Zq/IO) and population 
I (representative of the metallicity in the Milky Way 
disk; Z > Zq/10) stars c ontinue to form in parallel 
can also be expected (e.g. [Salvaterra. Ferrara fc Daval 
[2010| ). an d signatures of this may already have been 
detected ([Jimenez fc Haiman|[2006l) . Even more exotic 
galaxies may be envisioned if the dark matter of the 
Universe has the properties required for the formation of 
long-lived population II I stars fueled by WI MP annihila- 
tions in minihalos (e.g. [Spolvar et all[2008D . Many such 
"dark stars" may then end up within the first galaxies 
during their hierarchical assembly, and possibly imprint 
detectable signatures in their spectra ([Zackrisson et al.l 
[20T0t ). In this paper, however, we focus on pure popula- 
tion III galaxies. 

Here, we present Yggdrasi0, a new population synthe- 
sis code custom-designed for modelling the spectral en- 
ergy distributions (SEDs) of high-redshift galaxies. To 
reflect the significant variance in terms of stellar content 
that these objects may display, Yggdrasil is equipped to 
handle mixtures of conventional population I, II and HI 
stars (hereafter pop I, II and HI) as well as dark stars. 
This code also allows the treatment of nebular emission 
from the photoionized gas and extinction due to dust. A 
number of Yggdrasil model grids are publicly available 
from the lead author's homepagc0. 

The SEDs of Single SteUar Populations (SSPfEl), from 
various other population synthesis models can be used as 
input to Yggdrasil, which will then reweight the SSP time 
steps to accommodate arbitrary star formation histories. 
Currently, t he SSP models of IZa ,c krisson et_al. (200_lj), 
Starb urst9 9 ([Leitherer et n][1999t IVazquez fc Leithered 
[2005f ) and [MarastonI ( 20051 ) are implemented as options 
for pop I and II stars. For the duration of this paper, 
we will adopt Starburst99 SSPs ge nerated with Padova- 
AGB stellar evolutionary t racks ([Vazquez fc Leithereri 
[2005( 1 and the Kroupal ([2001[ ) universal stellar initial mass 
function (IMF) throughout the mass range 0.1-100 Mq 
for population II (assumed metallicity Z = 0.0004) and 
pop I (assumed metallicity Z = 0.020) galaxies. 

While both theoretical arguments and numerical sim- 
ulations give strong support to the notion that pop HI 
stars must have been more massive th an the pop II and I 
stars forming later on (for a review, see [Bromm fc LarsonI 
[2004f ). the exact IMF of pop HI stars remains unknown. 
It has been argued that the Universe may have produced 
two classes of pop HI stars - pop III.l stars which formed 
first, with characteristic masses of about ^ 100 Mq, 
and pop HI. 2 stars which formed somewhat later and 
had lower characteristic masses of ~ 10 Mq due to HD 
cooling promoted by the Lyman- Werner feedback pro- 
vided by the pop III .l stars (e.g. [Mackev et all [2003[ : 
[GreiffcBroI^ [20061 ). Naively, one would expect the 
pop III.l IMF to be appropriate for the stars forming 
in isolation (or in small numbers) within ~ 10^^^ Mq 



° named after a sacred tree in Norse mythology 
^ www.astro.su.se/~ez 

also known as single-age stellar populations or instantaneous- 
burst populations 
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Fig. 1. — The nebular contribution to the rest-frame SEDs of young, dust-free a) pop I (Z = 0.020) and b) pop III.l galaxies. Blue lines 
represent the purely stellar SEDs, whereas the red lines represent the total SEDs (including both stellar and nebular contributions). An 
instantaneous burst and an age of 1 Myr is assumed in both cases. The nebular component completely transforms the appearance of the 
SED in both panels, but the relative contribution of nebular emission becomes much higher for the pop III.l object (right panel) because 
of a higher fraction of hot, high-mass stars, which boost the Lyman continuum at A < 912A. The stellar SEDs have been rescaled to the 
same (arbitrary) flux level at 1300 A to facilitate a comparison between the panels. 



minihalos capable of H2 cooling at z « 20-40, whereas 
the pop III. 2 IMF may be more relevant for the first 
pop III galaxies forming in > 10'^~^ M(7) halos capable 
of HI cooling at 2 < 15 (e.g. iJohnson. Greif fc Bromrnl 
I2008| ). However, this pop III.l (~ 100 Mq) versus pop 
HI. 2 10 Mq) picture is looking increasingly uncer- 
tain, as the latest simulations indicate more fragmen- 
tation t han before , with lower pop HI masses a s a re- 
sult (e .g. IStacv et a l. 2008; Clark et al. 2011; Grei f et al.l 
l2011al) . It has recently also been recognized that su- 
personic streaming veloc ities of baryons with respect 
to the dark mat ter (e.g iTseliakhovich fc Hiratal 120111 : 
iMaio et all 120111 ) could have increased the turbulence 
within minihalos a nd assisted in re ducing the masses 
of pop III.l stars (|Greif et al.l 1201 lb). For simplicity, 
we have chosen to stick to the traditional pop III.l 
100 Mq) and HI. 2 10 Mq) convention through- 
out this paper. To cover all the bases, we adopt both 
of these top-heavy IMFs as viable options for our pri- 
mordial gala xies , but also consider pop HI galaxies with 
the same Kroupa; (2001) IMF as that used for pop II/I 
galaxies. 

For pop III.l galaxies, we will adopt the iSchaereil 
(pool stellar SSP with a power-law IMF {AN/dM cx 
M~") of slope a = 2.35 throughout the mass range 
50-500 Mq. For po p HI. 2 galaxies, we will adopt the 
iRaiter et al.l ([20 10b) TA model, which has a log-normal 
IMF with characteristic mass Mc = 10 Mq, dispersion 
a — 1 Mq but wings extending from 1-500 Mq. 

The contribution to the SED from photoionized 
gas is computed usi ng the procedure outlined in 
iZackrisso n et al.l ()2001[ ). In this machinery, the stellar 
population SED is, at every age step (once the star for- 
mation history has been taken into accou nt), used as in- 
put t o the photoionization code Cloudy (|Ferland et al.l 
[l99l . This resuhs in a self-consistent prediction for 
the nebular continuum and emission line fluxes which 
reflects the temporal changes in the ionizing radia- 
tion field. Other approaches in the literature include 
calculating the nebular continuum based on tabulated 
emissivities assuming a fixed electron temperature (e.g. 



iLeitherer et"an I1999L see IRaiter et all l2010bl for a de- 
scription of the shortcomings of this technique), or using 
emission-line ratios fixed by ern pirical calibrations (e.g. 
I Anders fc Fritze-v. Alvensleben|[2003 V While the latter 
technique may be useful for certain classes of objects, it 
will of be of little use for objects with properties that 
differ significantly from the empirical templates, since 
the emission-line ratios are expected to vary as a func- 
tion of age, IMF, metallicity and the physical conditions 
in the nebula. Throughout this paper, we will assume 
the nebula to be spherical, ionization-bounded and to 
have a constant, non-evolving hydrogen density n(H). 
These assumptions concerning the properties of the neb- 
ular gas, and their impact on our conclusions, are fur- 
ther discussed in Sect. |3l Additional Cloudy parameters 
include the gas filling factor /fin (which describes the 
porosity of the gas), the gas covering factor /cov (which 
regulates the amount of Lyman continuum leakage into 
intergalactic space) and the gas metallicity Z^^s- 

In the case of plane-parallel nebulae, there is a homol- 
ogy relation between photoionization models with the 
same ionization parameter U : 



U = 



Q(H) 



47rci?2n(H) ' 



(1) 



where (5(H) is the number of ionizing photons per unit 
time and R is the distance to the ionization source. The 
situation becomes more complicated for spherical nebu- 
lae (as assumed here), due to the dilution of the radiation 
as it propagates through the nebula (giving high U close 
to the centre and low U further out). Instead of specify- 
ing U directly, the inner radius of the cloud, Rm and the 
mass available for star formation, Mtot are instead used 
as input to Yggdrasil. The latter is defined as the gas 
mass converted into stars during a star formation episode 
of duration r: 



Mtot, = 



SFR(t) dt. 



(2) 



While the shape of the stellar SED is set by the stel- 
lar SSP, the age and the star formation history, Mtot 
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determines the overall scaling, and hence luminosity and 
(3(H). Please note that Mtot does not represent the mass 
of the resulting nebula, which may either be higher or 
lower than Mtot, depending on (5(H) and the gas param- 
eters (n(H), /fill, /f-nv and to smalle r extent ^gas)- 

As in the iZackrisson et al.l (|2001| ) model, the inner ra- 
dius of the ionizing cloud is set to: 

/ r \ 1/2 

i?in = 100 Rq i—j , (3) 

where L represents the bolometric luminosity of the 
model galaxy. This gives an effective mean ionization 
paramet er in the same range as observed in local HII 
regions C Kewlev fc Do^ [2001 . if fmi = 0.01 and 
n(H) = 100 cm~'^. Throughout this paper, we adopt 
Mtot — 10® Mq, to reflect the likely stellar population 
masses of the faintest galaxies detectable with JWST (see 
Sect. S]). However, we have verified that the resulting 
JWST colours remain practically unchanged for objects 
that are a up to a factor of ~ 10^ more massive. 

While dust extinction is not expected to be present 
during the very first star formation episode of pop HI 
galaxies, it may well become important after a few Myr, 
when pair-instability supernovae from 140-260 Mq stars 
or type II supernovae with M < 50 (mass lim- 
its valid in the absence of stellar rotation; i Heger et al.l 
[200 2') release their nucleosynthesis yields into the sur- 
roundings. In pop II/I galaxies, extinction is likely to 
be relevant at all ages. To allow for the treatment of 
dust extinction, Yggdrasil allows a choice of four diffcr- 
ent dust-c orrect ion recipes: the Milky Way, LMC, SMC 
(|Pe]||1992D and ICalzettI (|1997D attenuation models. In 
the latter case, corrections are applied separately to the 
nebular and stellar contributions to the SED, so that 
the nebular emission experiences a higher extinction. All 
dust extinction corrections are applied after the nebular 
SED has been generated. This is equivalent to assuming 
that the dust is located outside the HII region, and hence 
does not affect the Lyman continuum flux prior to gas 
absorption. This is reasonable assumption, at least for 
young pop HI galaxies experiencing a brief star formation 
episode, since no current models predict dust formation 
in the immediate vicinities of pop HI stars. 

3. THE NEBULAR PROPERTIES OF HIGH-REDSHIFT 
GALAXIES 

Nebular emission typically has a strong impact on the 
SEDs of young or star-forming galaxies, and the rela- 
tive contribution from nebular gas is expected to become 
stronger if the metallicity becomes very low or if th e IMF 
turns m ore top-heavy (Schacrcr 2002, 2003; Raiter ~e"t al.l 
I2010bf) . This is demonstrated in Fig. [TJ where we com- 
pare the rest-frame SEDs of 1 Myr old SSP models for a 
pop I galaxy (Z = 0.020. lK7mipall2Q0ll IMF) and a pop 
III.l galaxy. In both cases, the blue line represents the 
purely stellar SEDs whereas the red line describes the 
total SED, i.e. the SED including both stellar and neb- 
ular contributions. The nebular component completely 
transforms the appearance of the SED for both types of 
galaxies, but the relative contribution of nebular emis- 
sion becomes much stronger for the zero-metallicity, top- 
heavy IMF object (pop III.l). 

As we have argued in several previous papers 



(e ^Zackrisson et al.|[2008t iSchaerer fc de Barroi 120091 
I2010D . the flux boost due to photoionized gas can also 
have a notable impact on the broadband fluxes of high- 
redshift galaxies. In Fig. [2^, we demonstrate this by 
plotting the ratio of nebular to stellar flux, /neb //stars, in 
the JWST/NIRCam F444W fiher (at 4.44 ^m) as func- 
tion of redshift (within the range z — 0-20) for new- 
born (age 1 Myr), galaxies with SEDs representative of 
pop I (Z = 020), pop II (Z = 0.0004), pop HI with 
iKroupal (pOOl IMF, pop III.2 and pop III.l SEDs, as- 
suming an instantaneous burst. Nebular emission is seen 
to dominate the flux in this filter at all redshifts con- 
sidered for the pop HI models. As expected, the highest 
/ncb//stars are produced by the pop III.l IMF (cyan line), 
f ollowed b y the pop HI. 2 IMF (blue line) and pop HI with 
alKroupal (2001) IMF (green line). The pop II (red line) 
and pop I models (black line) produce significantly lower 
ionizing fluxes per unit mass. 

The bumps and wiggles along the lines in Fig. [5^ are 
produced when various emission lines and continuum fea- 
tures redshift in and out of the filter. The bump at z w 5- 
7 is due to the Ha line at 6563 A, and the extension of 
this bump to z « 9 in the case of the pop II galaxy 
is primarily due to the [OIII] line at 5007 A, which is 
strong at low to intermediate metallicities, weak in high- 
metallicity objects and completely missing from pop HI 
spectra due to lack of metals in these objects, as further 
discussed in Sect. [5j The overall increase in /neb //stars 
towards very low redshifts is due to the very different 
spectral slopes of the stellar and nebular continuum at 
rest- frame wavelengths above 10^ A (see Fig. [T]). 

In Fig. [^b, we plot the /ncb//stars ratio at z = 10 as 
a function of age for the same instantaneous-burst mod- 
els as in Fig. [2ti.. At an age of 1 Myr, going from a 
IKroupal (|200lD IMF pop I (black hue) to a pop II (red 
line) increases the /neb//stars ratio by a factor of « 2.5, 
and going from a pop II to a pop HI galaxy (green line) 
boosts /neb//stars by an additional factor of w 4. Shifting 
to more top-heavy IMFs (blue for pop HI. 2 and cyan for 
pop III.l) raises /neb//stars evcu further. These overall 
ratios between the models are approximately stable as 
a function of age, except for the pop III.l model, which 
evolves much faster than the others due to the lack of 
stars with masses below 50 Mq, fading from sight after 
just 3 Myr. I n the case of the pop HI. 2 and pop HI, 
IKroupal ()2001D IMF models, it takes « 12-16 Myr before 
direct star light starts to dominate over nebular emis- 
sion (thin horizontal dashed line), and « 40-50 Myr un- 
til nebular emission gives a negligible contribution (10%, 
as indicated by the thick horizontal dashed line) to the 
flux. 

Even though Lyman continuum leakage from high- 
redshift galaxies may well reduce the relative contribu- 
tion from nebular emission (as further discussed below), 
it is clear from Fig. [2] that it would take a very large 
amount of leakage to break the dominance of nebular 
light in the SEDs of young pop HI galaxies. For a 1 Myr 
pop III.l galaxy, /ncb//stars w 20 at z = 10, whereas the 
corres ponding values for pop HI. 2 and pop HI galaxies 
with a lKroupal (2001) IMF are /nob//stars ~ 15 and w 10 
respective ly. After 10 Myr, pop HI galaxies with a pop 
III.2 or a IKroupal pOOl) IMF have /ncb//stars ~ 5 and 
« 3, respectively. To bring nebular emission down to 
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Fig. 2. — The impact of nebular emission on pop I, II and III galaxies, a) The ratio of nebular to stellar light as a function of redshift 
for a 1 Myr galaxy in the JWST/NIRCam F444W filter at 4. 44 iim. The di fferent colours correspond to different assumptions about the 
metallicity and IMF: Pop I (black), pop II (red), pop III with a lKroupal l|2001i l IMF (green), pop III. 2 (blue) and pop III.l (cyan). The pop I 
and pop III.l models correspond to the SEDs depicted in Fig.[Tl The bumps and wiggles along the lines are produced when various emission 
lines and continuum redshift in and out of the filter. Nebular emission is seen to dominate the observed fluxes (logj^g /neb//stars > 0; thin 
horizontal dashed line) at essentially all redshifts for Pop III and Pop II, but the exact contribution depends on the metallicity and IMF. 
The relative flux contribution from the nebula is always more than an order of magnitude higher in the population III.l galaxy than in the 
population I galaxy, b) The evolution of /ncb//stars with age in the JWST/NIRCam F444W filter for instantaneous burst models. The 
line colours have the same meaning as in the left panel. The pop III.l predictions (cyan line) ends abruptly at w 3 Myr, since no stars 
in the 5O-5OOM0 range have life spans longer than this. In the case of the pop III. 2 and pop III, Kroupa (2001) IMF models, it takes 
^ 12—16 Myr before direct stars start to dominate over nebular emission (logj^Q /ncb//stars < 0; thin horizontal dashed line), and ^ 40-50 
Myr until nebular emission gives a negligible contribution (logj^Q /ncbZ/stars < —1; thick horizontal dashed line) to the flux. In the case of 
young (< 10 Myr) pop III galaxies, it would take a very large amount of Lyman continuum leakage (/osc > 0.95) to bring nebular emission 
down to a level where it no longer affects the broadband fluxes in any significant way. 



a level where it no longer affects the broadband fluxes 
of < 10 Myr old pop III objects in any significant way 
(/ncb//stars 0.1), the ionizing escape fraction would 
have to be /osc ^ 0.95. 

It is doubtful whether pop III galaxies can be identified 
based on colour criteria for much longer than ^ 10*" yr. 
The lifetimes of massive pop III stars are very short (w 2- 
3 X 10^ yr years for « 50-500 Mg stars; [Schaercr 2002), 
and as soon as the first supernovae explode, metals will 
be dispersed into the surrounding medium. While the 
ejecta from the first pop III supernovae may take ~ 10® 
yr to cool sufhciently to be us ed in the subsequ ent forma- 
tion of pop II and pop I stars CG reif et al.ll2007l ) , the time 
it takes before metal emission lines start to emerge from 
the surrounding nebula may be much shorter, perhaps no 
more than a few Myr. At that point, the broadband sig- 
natures of metal- free nebulae discussed in Sect. [5] would 
be jeopardized. If feedback from the first supernovae 
clears the galaxy of photoionized gas, the SED would be 
dominated by direct star light from pop III stars until 
pop II/I star formation sets in (after ^ 10® yr), but even 
then, the unique spectral characteristics of pop III stars 
are only retained for ^ 10^ yr (see Sect. [6]). 

3.1. Type A, B and C 

While Fig. 11] and Fig. [5] would suggest that the SEDs of 
pop III galaxies are almost completely nebular in nature, 
there are a number of mechanisms that could diminish 
the imp act of nebular emissio n. As shown in the simula- 
tions bv lJohnson et al.l (|2009f) . stellar feedback can drive 
the photoionized gas present in these systems outward, 
eventually pushing it beyond the virial radius of the dark 
matter halo. Once the HII region breaks out of its host 
halo, the gas density will be extremely low, resulting 



in a huge intergalactic nebula with very long recombi- 
nation timescales. This nebula, while very important 
for cosmic reionization, is probably of little importance 
for JWST observations of individual galaxies. When at- 
tempting aperture photometry of galaxies surrounded by 
such large, low-surface brightness nebulae, the nebular 
contribution (appearing in both the object and sky an- 
nulus) may be almost completely subtracted away, effec- 
tively leading to the detection of a SED dominated by the 
more concentrated stellar component. Whether or how 
quickly the gas gets expelled depends on the star forma- 
tion efhciency within a given halo, and also the form of 
the IMF. Having a very low ratio of stellar to dark mat- 
ter mass, and/or an IMF without too many high-mass 
stars, would likely result in a situation in which the HII 
region stays completely confined within the CDM halo 
(for the limiting case of just one pop III star in h alos 
with different masses, see e.g. iKitavama et al.ll2004[) . 

In Fig. [3l we have schematically defined three classes 
of objects in different stages of gas expulsion. For type 
A, the ionization-bounded HII region (grey region) re- 
mains compact and confined well inside the virial radius 
(dashed line). This implies a Lyman-continuum escape 
fraction /osc ~ and a maximal contribution from neb- 
ular emission to the overall SED of the galaxy. In the 
case of type B, stellar feedback has pushed the HII- 
region partly outside the virial radius {'^ 1 kpc for a 
10® Mq halo at z « 7-13), giving rise to partial leak- 
age of ionizing radiation into the intergalactic medium: 
< /esc < 1- For type C galaxies, the ionizing gas 
has been expelled outside the dark matter halo, result- 
ing in an SED with minimal nebular contribution and a 
Lyman-continuum escape fraction /osc ~ 1- While this 
simple picture admittedly neglects the anisotropic out- 
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Fig. 3. — Schematic illustration of the dynamical evolution of the ionized gas surrounding the first galaxies, a) For type A galaxies, 
the HII region (dark region) is compact, ionization/radiation bounded and confined well within the virial radius (dashed fine) of its dark 
matter halo. This implies a very low escape fraction for ionizing radiation (/esc ~ 0) and a maximal contribution of nebular emission to 
the overall SED of the galaxy, b) In the case of type B, part of the HII region has been pushed outside the virial radius of the halo, with 
partial escape of ionizing radiation as a result (0 < /esc < !)• At this point, the SED may still contain signatures of nebular emission, but 
at a lower level than that provided for type A. c) In type C objects, the HII region has been completely pushed outside the virial radius, 
implying a massive escape of ionizing photons into the intergalactic medium (/esc ~ !)• At this point, the nebula will be huge, have very 
low density, long recombination timescales and probably very little impact on the observations of individual galaxies. Hence, the SED is 
expected to be almost completely stellar, with essentially no contribution from nebular emission. 



flow of gas and the irregularities within th e photoionized 
medi um evident from actual simulation ([Johnson et alj 
|2009() . it still captures the salient points when it comes 
to modelling the nebular contribution to the galaxy SED. 

In this paper, we will only treat galaxies that belong to 
categories A (full nebular contribution) and C (no neb- 
ular contribution). The intermediate type B would be 
far more challenging to model in detail - the gas density 
profile is highly time-dependent at this stage, the long 
recombination timescales introduce a time lag between 
the evolution of the stellar and nebular SEDs, and com- 
monly used sky subtraction strategies may interfere with 
the predicted contribution from the the nebula to the ob- 
served SED. While the relative flux contribution of the 
nebula - and therefore the overall luminosity - does de- 
pend on these effects, the colours of type A and type B 
galaxies are nonetheless likely to be similar as long as the 
system is young and nebular emission is dominant. The 
youth criterion comes in because the mismatch between 
the momentary stellar SED and the prior ionizing field 
to which the nebula responds becomes more severe once 
the high-mass stars of the stellar population have died 
and many intermediate-mass stars have evolved off the 
main sequence. 

The nebular contribution from type A galaxies is mod- 
elled assuming a constant-density HII region with hy- 
drogen density 7t,(H) — 100 cm^'^ and filling factor 
/fill = 0.01. In the case of pop III galaxies, the gaseous 
metallicity is set to Zg^s = 10~^ (whereas Zgtars = 0), 
whereas pop II {Z = 0.0004) and I {Z = 0.020) galax- 
ies have Zgas = -^stars- Scalcd solar abundances are 
used in all cases. The default value for the gas cover- 
ing factor of type A galaxies is /cov = 1 (implying a 
Lyman continuum escape fraction f^sc — 0). However, 
in the case of a compact HII region, anisotropic feed- 
back, supernova chimneys and irregularities in the gas 
density may result in large holes in the nebula, through 
which ionizing radiation can escape into the intergalac- 
tic medium. Simulations predict that the amount of Ly- 
man continuum escape ma y be a function of dark halo 
mass (Gncdin ct al. 2008; Razoumov fc Sommer-LarsenI 
120101: lYaiima et al...20111 but the likely escape fraction 
and its exact mass dependence remains controversial. In 



TABLE 1 

The mAB magnitiude limits adopted for, lOo- detections 

AFTER texp = 100 H EXPOSURES, BASED ON THE NIRCAM AND 

MIRI sensitivities listed on the JWST homepage''. 



Instrument 


Filter 


A (/im) 




NlROam 


F070W 


0.70 


29.9 




F090W 


0.90 


30.4 




F115W 


1.15 


30.6 




F150W 


1.50 


30.7 




F200W 


2.00 


30.9 




F277W 


2.77 


30.6 




F356W 


3.56 


30.6 




F444W 


4.44 


30.1 


MIRI 


F560W 


5.60 


27.6 




F770W 


7.70 


27.3 




FIOOOW 


10.0 


26.3 




F1130W 


11.3 


25.3 




F1280W 


12.8 


25.5 




F1500W 


15.0 


25.2 




F1800W 


18.0 


24.3 




F2100W 


21.0 


23.5 




F2550W 


25.5 


22.3 



http: / /www. stsci.edu/jwst/instruments/ 



our model, the possibility of Lyman continuum escape 
through holes is treated by allowing the gas covering fac- 
tor /cov to take on values < /cov < 1- While having a 
very low but non-zero /cov does not capture all the com- 
plexities of a type B galaxy, it still gives some indication 
of what to expect from such an object. 

Type C objects are treated by setting /cov = 0, which 
implies a Lyman-continuum escape fraction of unity and 
an SED composed of stars only. Since no nebula is pro- 
duced, a photoionization code like Cloudy is not required 
to model this case. 



4. THE DETECTION LIMITS OF HIGH-REDSHIFT 
GALAXIES IN JWST BROADBAND SURVEYS 

In Fig. m we use Yggdrasil to predict the popula- 
tion masses of the faintest star-forming galaxies de- 
tectable through JWST broadband imaging when tak- 
ing extremely long exposures of a small patch of the sky, 
i.e. a so-called Ultra Deep Field (UDF) . The mass limits 
presented are based on the requirement that galaxies are 
detected at lOcr in at least one JWST broadband filter 
(spectral resolution R = i) after a 100 h exposure (per 
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filter). While the detailed specifications of JWST UDFs 
have not yet been fixed, 100 h represents a reasonable es- 
timate of the longest exposure times that are likely to be 
relevant. The resulting mass limits can easily be rescaled 
to other exposure times toxp using: 



Mr, 



oc 



in 



100 h 



1/2 



(4) 



Our derived mass limits take all of the 17 JWST broad- 
band [R — 4) filters into account, and are based on the 
actual JWST transmission curves of these filter^^ and 
the estimated AB-magnitude limits listed in Table. [TJ 



These lOcr, to 



100 h detection limits rriAB, looh can 



be rescaled to other exposure times using: 



mAB = 2.5 log 



10 



100 h 



1/2- 



"lAB, lOOh- 



(5) 



To avoid predictions hinging on the highly uncertain 
luminosity of the Lyg emission line at h igh redshifts (e.g. 
iFontana et all 12010: IHaves et al.ll20Tll) . we assume the 
Lyof escape fraction to be zero - i.e. the Lya line does 
not contribute to the predicted fluxes at all. When com- 
puting the JWST broadband fluxes, we furthermore set 
all SED fluxes shortward of Lya to zero for model galax- 
ies at z > 6, to reflect the high level of absorption in the 
neutral intergalactic medium at these epochs. 

The mass limits in Fig. 2] are presented for pop I, pop 
II, pop III galaxies of both type A (maximal contribu- 
tion from the nebula) and C (no contribution from the 
nebula) with constant star formation rates as a function 
of age (up to 30 Myr) and redshift (up to z = 20). As 
discussed in Sect. [3l galaxies starting from zero metal- 
licity are not likely to retain their pop III characteristics 
for much longer than this. While pop I and pop II galax- 
ies may well have higher ages at these epochs, such ob- 
jects would need to be even more massive than suggested 
by the limits at 30 Myr in Fig. |4] to be detectable with 
JWST. Hence, Fig. |4] places hard limits on what can re- 
alistically be observed with this telescope through ultra 
deep imaging in the absence of gravitational lensing. 

The mass used to define these limits is the mass of 
gas converted into stars since the beginning of a star 
formation episode: 



M{t) 



SFR(t') At'. 



(6) 



This is equivalent to the population mass often dis- 
cussed in the context of simulations, where the overall 
gas mass is multiplied by the star formation efficiency of 
the first starburst episode to compute the gas converted 
into stars. The mass in luminous stars (stars that have 
not yet exploded as supernovae or turned into compact 
remnants) at these ages can be considerably lower at ages 
> 3 Myr (especially in the case of top-heavy IMFs) , since 
many of the stars forming at t = yr have then already 
faded away. For instance, the mass in luminous stars in 
a pop III.l galaxy with constant star formation is lower 

The one exception is the NIRCam/F090W filter profile, which 
was not available to us at the time of writing. In this case, the 
results are instead based on a top-hat transmission curve. 



than the plotted limits by a factor of w 5 after 10 Myr, 
and by a factor of 14 at 30 Myr. On the other hand, 
the mass in luminous stars within IKroupal (|200lD IMF 
populations is reduced by only « 20 % (pop III) and 
« 10 % (pop II and I) after 30 Myr of constant star 
formation. In the case of star formation episodes with 
finite duration r, M{t) equals Mtot defined in eq.® if 
evaluated at a time after the end of the starburst (t > t). 

As seen in Fig.jH M ~ 10^ Mq star-forming pop I and 
pop II galaxies (black and red lines) can be detected at 
z « 10, whereas pop III.l galaxies (cyan lines) can be 
detected even if the mass is as low as M 10^ Mq. Pop 
III. 2 galaxies (blue lines) lie between pop I/II and pop 
III.l galaxies in this dia g ram. On the other hand, pop III 
galaxies with a IKroupal (|200lD IMF (green lines) tend to 
require a mass at least as high as pop I/II galaxies (and 
sometimes significantly higher) to be detected, because of 
their lower stellar mass-to- light ratios in these passbands. 

For galaxies with a given IMF and metallicity. Type 
A objects (maximal nebular flux contribution) are al- 
ways brighter than type C objects (no nebular flux con- 
tribution) and can therefore be detected at lower masses. 
However, the difference is much smaller for pop I and II 
galaxies (less than a factor of 2 at 2 = 10) than for pop 
His, due to the higher relative ionizing fluxes of the lat- 
ter objects. For instance, nebular emission allows young 
(« 3 Myr old) pop lll.l galaxies at z = 10 to be detected 
at masses a factor of « 4 lower than in the case where 
all the ionizing radiation is escaping into the intergalactic 
medium. The mass detection limits presented for pop HI 
and pop II galaxies in Fig. [3] are similar (within a factor 
of 2, once differences in exposure times are considered) 
to those derived by Pawlik et al. (2011), who considered 
the JWST detectability of lines and continuum radiation 
separately. 

The detection limits in Fig. |3] indicate the masses of 
galaxies that JWST would be able to detect at a single 
wavelength, but to derive the properties of the objects 
observed, photometry in several filters is always required. 
Since the JWST sensitivity varies greatly as a function of 
wavelength, the detection limits are sensitive to the fil- 
ters used. In Fig.jSl we plot the lowest stellar population 
masses of 3, 10 and 30 Myr old type A galaxies (maxi- 
mal nebular emission) detectable in a JWST UDF (lOcr 
detections after 100 h exposures) as a function of central 
filter wavelength at z — 6, 8, 10 and 12. The mass lim- 
its stay approximately constant within the wavelength 
range of the JWST/NIRCam instrument (operating at 
0.7-4.4 /im), but then rise by an order of magnitude as 
one switches to JWST/MIRI observations at > 5.5 fj,m 
and continue to deteriorate almost monotonically all the 
way up to the 25 /im limit of MIRI. Hence, detections 
of high-redshift galaxies in the MIRI filters will be very 
expensive in terms of observing time, compared to those 
in NIRCam filters. On average, the relative contribu- 
tion from nebular continuum grows with increasing rest- 
frame wavelength, and this compensates for some of the 
sensitivity loss in the MIRI bands. For galaxies without 
nebular emission (type C), the decline in mass sensitiv- 
ity with increasing JWST passband wavelength is even 
more dramatic than plotted here. 

5. THE SPECTRAL SIGNATURES OF TYPE A, POP III 
GALAXIES 
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Fig. 4. — Predicted mass detection limits in JWST ultra-deep fields. The panels show the lowest masses of burst-like stellar populations 
that JWST can detect through broadband imaging at lOcr after a 100 h exposure, as a function of redshift. The l ine colours rep resent 
different combinations of population metallicities and IMFs: Pop III.l (cyan line), pop III. 2 (blue), pop III with a IKroupal 1 120011 ) IMF 
(green), pop II (red line) and pop I (black line). The left column contains the predictions for type A galaxies (maximal nebular contribution) 
and the right column the corresponding predictions for type C galaxies (no nebular contribution). Each row corresponds to a different 
starburst age (3, 10 and 30 Myr from top to bottom) for a population undergoing a star formation episode with constant star formation 
rate. 



The JWST fluxes and colours of pop III galaxies 



9 



z=6 z=8 




5 10 15 20 25 5 10 15 20 25 

Wavlength (|im) Wavlength ()xm) 



Fig. 5. — Predicted mass detection limits for young, type A galaxies in JWST ultra-deep fields as a function of filter wavelength. The 
panels show the lowest stellar population masses of burst-like, 10 Myr old, type A (maximal nebular contribution) galaxies that JWST 
can detect through broadband imaging at lOcr after 100 h exposures per filter, as a function of central filter wavelength. The line colours 
have the same meaning as in Fig. U Each panel corresponds to a different redshift (z = 6, 8, 10, 12). No results are included for filters at 
< 1.15 fira, since IGM absorption is likely preventing detection in these filters at some of these redshifts. The sharp increase in the lowest 
detectable mass for the FISOW filter (at 1.50 iim) at z = 12 is also due to this effect. The dashed vertical line in each panel delimits the 
wavelength ranges of the NIRCam to MIRI instruments. 
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As discussed in Sect. [2 all searches for pop III 
galaxies with largely nebular SEDs (our type A) 
have so far focused on the fluxes of the Lya line, 
the Hell A1640, A4686 lines, or the Lyman 'bump ' 
(e.g. [T umlinson fc ShuU' 2000; Tumlinson ct al. 2001; 
Oh et al. 20qi:.Malhotra fc Rhoads 2002; Schacrcr 2002, 
20031 : llnouel I2O1O0 . However, neither of these spectral 
features is likely to produce a robust and detectable sig- 
nature in broadband filters for objects at z > 6. Hell is 
too weak, the Lyman 'bump' is absorbed by the inter- 
galactic medium and the escape fraction of Lya photons 
is notoriously difficult to predict. 

Here, we propose a different strategy, based on the 
fact that such galaxies are likely to have strong hydro- 
gen and helium emission lines, but - due to the lack 
of metals - no lines due to heavier elements. Cer- 
tain rest frame UV and optical metal emission lines, 
most notably [OII]A3727, [OIII]A5007 and [SIII]A9069, 
are sufficiently st rong to influence the broadband fluxes 
of galaxies (e.g. lAnders fc Fritze-v. AlvenslebenI 120031 : 
iZackrisson et al.l 120081: ISchaerer fc de Barrosll20'09l) . and 
removing these lines will therefore al ter the broadband 
colours of type A objects. Recently, llnouel (|2011b[ ) ar- 
gued that the lack of [OIIIJA5007 emission would be one 
of the best ways to single out pop HI galaxies. Our anal- 
ysis agrees with this conclusion. 

5.1. NIRCam colour criteria 

'Inoue" (2011b) suggests that imaging in the two 
JWST/NIRCam fihers F277W and F444W could be use- 
ful to single out pop HI galaxies at z « 8, since the 
lack of [OIII] emission in the F444W filter would ren- 
der the colours of such objects very blue. However, 
the m277 — TO444 < criterion proposed by Inoue for 
pop HI objects is not unique, as demonstrated using 
instantaneous-burst models in Fig. |6] ~ both pop II and 
pop I galaxies of type A can at certain ages venture blue- 
ward of this limit. In the case of our pop II model 
(red line), this happens because of a short-lived evo- 
lutionary phase at ages of ~ 10 — 20 Myr, where the 
ionizing flux has dropped sufficiently to make the neb- 
ular contribution irrelevant, yet the UV continuum of 
the stars remains very blue. In the case of our pop I 
model {Z = 0.020; black solid line), this happens be- 
cause very young (< 5 Myr) populations with high metal- 
licitics also display very we ak [OIII] emission lines (e.g. 
Hanagia 2003; Na gao et al.l [20061. In an instantaneous 
burst scenario, the timescales over which which Pop I/II 
systems display these very blue colours are shorter than 
the corresponding timescales for pop HI systems, which 
may limit the risk of pop I/II interlopers in surveys for 
pop HI galaxies. On the other hand, the typical star 
formation histories may be different in the two cases, 
and adopting a more extended star formation history 
for pop I/II could make these objects display pop Ill- 
like colours for a longer period of time. Dust extinc- 
tion, which is likely to be associated with star formation 
in high-metallicity environments, may of course render 
the colours of young pop I/II systems redder and possi- 
bly prevent metal-enriched galaxies from entering colour- 
selected pop HI samples. 

While the details of the IMF determines the longevity 
of pop HI galaxies, the 771277 — TO444 colours predicted for 
type A, pop HI objects of age < 10 Myr are not sensitive 
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Fig. 6. — The m277 — 01444 colours of type A, instantaneous-burst 
galaxies at 2 = 8. The line colours indicate different metallicity and 
I MF comb inations: Pop III.l (cyan), Pop III. 2 (blue), Pop III with 
a [Kroupa (200T) IMF (green), Pop II (red) and Pop I (black). A 
Eroupa (2001) IMF has been assumed for all of the Z = 0.0004- 
0.020 populations, ilnoua (I2011ti ) suggests that m277 — 771444 < 
(shaded region) should serve as a colour criterion for pop III 
galaxies, but as seen, both Pop II (red line) and Pop I (black solid 
line) objects can meet this colour criterion during certain phases 
of their evolution. While the details of the IMF determines the 
longevity of pop III galaxies, the m277 — m444 colours predicted 
for objects of age < 10 Myr are not sensitive to the exact pop III 
IMF, as illustrated by the overlapping cyan, blue and green lines. 
The age axis runs up to the age of the Universe at this redshift 
6.5 X 10** yr). 

to the exact pop III IMF, as illustrated by the largely 
overlapping cyan, blue and green lines in Fig. |6l 

5.2. Combined NIRCam and MIRI colour criteria 

In Fig. [3 we demonstrate that that cleaner selec- 
tion criteria for type A, pop HI galaxy candidates 
at z ss 8 can in principle be derived by combining 
JWST/NIRCam data in the F444W fiher with obser- 
vations in two JWST/MIRI fihers: F560W and F770W. 
The lack of [OIII]A5007 in F444W and the presence of 
Ha in F560W makes young pop HI galaxies redder in 
m444 - (Fig. [Tla,) , yet bluer in mseo - ^770 (Fig. [7b) 
than other types of galaxies. Due to the [OIII] degener- 
acy between pop III and pop I galaxies, both can attain 
similar 771444 — "7560 colours. However, since there is no 
similar degeneracy in 777560— 777770, pop HI objects appear 
in a corner of the 777559 ~ "^-770 vs. 777444 — 777560 diagram 
(Fig. [8]) where no other galaxies should appear. Pop HI 
objects remain in the upper left corner (grey region) up 
to w 15 Myr after formation in the case of an instan- 
taneous burst, and longer in the case of more extended 
star formation histories (in principle up to ~ 100 Myr in 
the case of a constant SFR, although such high ages are 
unlikely to apply, as discussed in Sect. [S]). The colours 
of pop HI galaxies with different IMFs are almost iden- 
tical at young ages. We have verified that these criteria 
hold for objects up to an age of ~ 15 Myr as long as the 
Lyman-continuum escape fraction is /osc ^ 0-5j irrespec- 
tive of the pop HI IMF (but higher /osc are allowed in the 
case of younger ages or more top-heavy IMFs) . Fig. |8]also 
illustrates the effect of dust reddening on these colour cri- 
teria. Even though substantial dust reddening (if at all 
possible in pop HI galaxies) could conceal pop HI galax- 
ies by making them appear as newborn, high-metallicity 
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objects, there is no risk of false positives - i.e. no other 
type of galaxies in our model grid can spuriously appear 
as pop III objects at this redshift. 

5.3. The redshift evolution of type A galaxy colours 

One major obstacle for any colour criteria applied to 
SEDs dominated by nebular emission, is that they evolve 
very quickly as a function of redshift. In Fig. |9] we 
plot the redshift dependence of the 771277 ~ 771444 and 
77^560 — '71770 colours for young (1 Myr old) pop III, II 
and I galaxies of type A. As seen, both colours evolve 
dramatically as a function of redshift. By contrast, the 
colours of equally young type C galaxies (i.e. with SEDs 
dominated by direct starlight; dashed lines) hardly evolve 
at all. Even though the 771277 " 771444 < criterion picks 
up on pop III galaxies at z w 6.7-8.5, it also selects 
young pop I and II galaxies of both type A and C at 
z < 5 and at z > 6.5. The risk of confusion can be lim- 
ited, but not completely removed, by adding photometry 
in other NIRCam filters and applying drop-out criteria. 
For instance, additional observations in the F090W and 
F115W filters should effectively limit the redshift range 
of the sample to z « 6.5-9.0 since objects in this redshift 
range are expected to appear as drop-outs in F090W, but 
not in F115W, due to IGM absorption shortward of Lya. 

A MIRI-based pop III criterion like TO560—TO770 < —1.3 
(shaded region in Fig. [T]) is not likely to generate nearly 
as many interlopers. This colour is useful as a pop III 
diagnostic over the limited range z fa 7-8 (shaded region 
in Fig. [S]) only, but within that range, pop III galaxies 
produce bluer colours than all other kinds oi Z < 0.020, 
z > I galaxies in our model grid. Even though we have 
here only plotted the colours of very young galaxies (1 
Myr old) , this holds for all interloper ages (up to the age 
of the Universe at each redshift considered) . While there 
is some room for confusion with z < 1 galaxies - which 
due to dust emission in the F560W filter (not taken into 
account in our models) may potentially produce similar 
colours in this diagram - such objects are not likely to 
display apparent magnitudes in the same range as high- 
redshift pop III galaxies. Hence, objects that end up in 
the upper left corner of Fig. [51 and have apparent magni- 
tudes in the range expected for high-redshift galaxies, are 
likely to be z « 7-8 pop III galaxies even in the absence 
of additional redshift constraints. 

5.4. Sensitivity issues 

While this scheme would seem to give a cleaner selec- 
tion o f pop II I galaxy candidates than that proposed by 
llnoud (|2011b[ ). it suffers from one obvious drawback: the 
lower sensitivity of MIRI implies a minimum mass for 
detection that is an order of magnitude higher than in 
the case where only NIRCam filters are used (see Fig. [5]). 
Even in the case of 100 h exposures per filter, the mass 
converted into stars would need to be on the order of 
~ 10^ Mq to allow detection in the F560W and F770W 
filters. The maximum mass that stellar populations con- 
sisting entirely of pop III stars can reach is unknown, but 
current simulations suggest that unenriched halos are un- 
likely to a ttain total masses in excess of M ~ 10^ Mq 
at z > 7 (jTrenti et all [20091 ). To produce - 10^ Mq 
worth of pop III stars, essentially all of the baryons in 
a M ~ 10® Mq object would need to be converted into 



stars (somehow evading negative feedback effects) in a 
limited amount of time (up to ~ lO'' yr), which seems 
highly unrealistic. By hunting for pop III galaxies be- 
hind lensing clusters with mag nification /i « 100 (MACS 
J0717.5-f 3745; e.g. iZitrin et ahj I200I iZackrisson eTal] 
l2010f) . the required stellar population mass can in prin- 
ciple be lowered to ^ 10^ Mq (implying a star forma- 
tion efficiency of e ~ 10"^). However, since gravitational 
lensing reduces the volume probed within a certain red- 
shift interval by a factor equal to the magnification it- 
self, the number density of pop HI galaxies would need 
to be very high to make this strategy successful. De- 
tailed lensing calculations would be required to assess 
the prospects of endeavours of this type (Zackrisson et 
al., in prep.). The alternative would be to hunt for pop 
HI galaxies using a coarse set co lour criteria based on 
NIRCam-only data (|Inouell2011b[ ), accept the possibility 
of pop I/II interlopers and then do spectroscopic follow- 
up with JWST/NIRSpec for all candidates. The expo- 
sure times required to detect a z w 8 pop HI object of 
a given mass in F560W or F770W with MIRI is actu- 
ally comparable to that required to reach the continuum 
level at w 4.5 /im and spectroscopically confirm the ab- 
sence of a [OIIIJA5007 line with NIRSpec. Hence, unless 
the number of candidate pop HI galaxies turn out to be 
very large (w 100; limited by the number of NIRSpec mi- 
croshutters) within a given field (3' x 3') at the relevant 
magnitudes, the NIRCam + NIRSpec strategy may be 
more economical in terms of telescope time. 

6. THE SPECTRAL SIGNATURES OF TYPE C, POP III 
GALAXIES 

If pop HI galaxies can experience sufficient Lyman con- 
tinuum leakage to render the nebular contributions to 
their SEDs negligible, such type C objects would stand 
out in JWST surveys because of their very blue rest- 
frame UV continuum slopes. This is demonstrated in 
Fig. [TUl where we compare the SEDs of 1 Myr old, 
instantaneous-burst pop HI, II and I galaxies at rest- 
frame wavelengths from 1250 to 4600 (the range ac- 
cessible for the NIRCam instrument for a z = 10 tar- 
get). The stellar continuum is seen to become progres- 
sively steeper when going towards lower metallicities and 
more top-h eavy IMFs - a fact previously noted by e.g. 
iBromm et al. (2001), Schacrcr (2003), Tumlinson et al.l 
igOOS), 'Schacrcr & Pello (2005),' Bouwens ct al. (2010*), 
iSniguchi ct al. (2010) and Raiter ct al. (2010b). How- 
ever, the difference between the pop HI. 2 and pop III.l 
galaxies appears to be undetectably small. As seen in 
Fig. [TJ inclusion of nebular emission would make the ob- 
served continua substantially redder, thus potentially al- 
lowing type C galaxies to be identified by their very blue 
broadband colours. The filters would, however, need to 
be carefully selected for the redshift range of interest to 
prevent small amounts of lingering nebular emission from 
interfering with this spectral signature. 

Basically, a photometric measurement of the rest-frame 
UV continuum should ideally use filters with as large 
wavelength separation as possible. However, since IGM 
absorption shortward of Lya likely precludes the use of 
certain NIRCam filters at z > 6, we have here selected 
771200 ~ 771444, as a trade-off between redshift range and 
UV slope sensitivity. In Fig. [11] we demonstrate that 
777200 — W444 should be able do distinguish type C, pop HI 
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Fig. 7. — The — mseo (a) and msgo — ^770 (b) colours of type A, instantaneous-burst galaxies at 2 = 8. The line colours represent 

different metallicity a nd IMF combi nations: Pop III.l (cyan), Pop III. 2 (blue), Pop III with a Kroupa (2001) IMF (green), Pop II (red) 
and Pop I (black). A lKroupal l|20011 ) IMF has been assumed for the two latter populations. In 01444 — "iseOi young (< 15 Myr old) pop 
III galaxies appear redder than galaxies of almost all other ages and metallicities due to the absence of strong [OIII] emission in the m444 
filter, yet strong H« emission in rrtsgo- The only notable exceptions are newborn (< 3 Myr old) pop I galaxies, which also display very 
red colours. In mseo — "1770, pop III galaxies appear bluer than all other objects at this redshift due to their high H« equivalent widths. 
While the details of the IMF determines the longevity of pop III galaxies, the colours predicted during their lifetimes are not sensitive to 
the exact pop III IMF, as illustrated by the largely overlapping cyan, blue and green lines. Shaded regions indicate the regions where pop 
III galaxies would bo identifiable. In both panels, the age axis runs up to the age of the Universe at this redshift (» 6.5 X 10^ yr). The 
error bars on the observed colours would be Ri 0.15 mag in the case of lOcr detections in the relevant filters. 

TABLE 2 

UV SLOPE-*- 13 FOR INSTANTANEOUS-BURST MODELS 




Motallicity/IMF 


Type2 
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Pop III.l 
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Fig. 8. — The signatures of type A, pop III galaxies in mseo — 
m77o vs. 771444 — "1560 diagram at z = 8. The lines correspond to 
a limited set of the instantaneous-burst models presented in Fig. [7] 
(to avoid cluttering): Pop III. 2 (blue line), Pop II (red line) and 
Pop I (black line) . Markers along the tracks indicate ages of 1 Myr 
(circle), 3 Myr (square) and 10 Myr (triangle). The tracks then 
continue up to the age of the Universe at this redshift (Ri 6.5 X 10* 
yr) . Arrows indicate how the colours of a 3 Myr old pop II galaxy 
would b e shifte d in the case of LMC extinction (short arrow) and 
ICalzettil 119971) extinction (long arrow) for rest frame E{B — V) = 
0.25 mag. In the case of Calzetti extinction, this refers to the 
extinction affecting by the stellar component SED, whereas that 
affecting the nebular component is higher by a factor 2.3. The 
reddening vectors are similiar (but not identical) for other ages and 
metallcities. In the absence of extinction, young pop III galaxies 
occupy a unique corner (shaded region) of this diagram. While 
extinction (if at all effective in pop III objects) could potentially 
shift pop III galaxies into regions occupied by young pop I objects, 
dust cannot make other types of galaxies mimic the intrinsic colours 
of pop III objects. This diagnostic diagram looks approximately 
the same throughout the redhsift range 2 = 7—8, as discussed in 
the main text. The error bars on the observed colours would be 

0.15 mag in the case of lOcr detections in the relevant filters. 



^ Derived using /3 — 4.29(Ji25 — ^^l6o) — 2 at 2: — 7 lIBouwens et all 
[2O70t) 

^ Type A: /esc — (maximal nebular emission) 
Type B: /esc = 1 (purely stellar SED). 

galaxies from all kinds of pop I and II galaxies through- 
out the redshift range z « 7.5-13. 

In Fig. [TTk . we demonstrate that any type C galaxies 
at z = 8 (solid lines) that display m2oo ~ 771444 < —1.2 is 
bound to be a pop III, type C system. Nebular emission 
in type A galaxies would only drive the colour redward. 
For instance, all newborn (« 1 Myr old) type A galaxies 
start off with m2oo — TO444 > —0.4 mag at this redshift. 
The (t(to200 ~ '71444) ~ 0.15 mag photometric error re- 
sulting from demanding a lOtr detection (the detection 
threshold assumed in Fig. 2]) in both filters should be suf- 
ficent to separate pop III.l and pop III. 2 galaxies from 
pop II and I, but woul d make i t difficult to distinguish 
pop III galaxies with a iKroupal (^2001) IMF from pop I 
or pop II galaxi es. The 771277 — 7)7444 colour suggested 
bv llnoug (|2011bf) as a diagnostic is not as suitable, since 
the predicted difference between pop III and pop II/I 
galaxies is smaller. 

In the instantanous-burst models used in Fig. [TTk . pop 
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Fig. 9. — The redshift evolution of the a) m277 — "1444 and b) msgo — ^1770 colours of 1 Myr old galaxies. The solid lines represent 
predictions for type A galaxies: Po p III. 2 (bl ue solid), Pop II (red solid) and Pop I (black solid). The dashed lines represent the corresponding 
predictions for type C models. A IKroupal 12001 ) IMF is assumed for both pop I and pop II objects. Since all pop III objects display 
very similar colours at this age, only the pop III. 2 predictions are included to avoid cluttering. Galaxies with SEDs dominated by nebular 
emission (type A; solid lines) are seen to display a far more dramatic redshift evolution than those with SEDs dominat ed by direct star 
light (type C; dashed lines). In the left panel, the 771277 — 711444 < criterion (thin horizontal dashed line) suggested bv llnou^ Il20ig 'l for 
pop III galaxies at 2 ?a 8 is seen to also pick up on many other types o f young type C galaxies over a wide range of redshifts. Here, the 
shaded region indicates the redshift interval 2 6.7-8.5, over which the IInou3 lj2010l 'l criterion holds for our pop III model galaxies. The 
right panel demonstrates that Pop III, type A galaxies at 2 7-8 (shaded region) should appear at m^ao — 1)1770 < —1.3 (thin horizontal 
dashed line), which is bluer than all other galaxy model in this diagram, and hence provides a cleaner selection of pop III candidates. 

for about 6 Myr, and Pop III. 2 and pop III.l would meet 
the same colour criteria for as long as star formation re- 
mains active. Of course, a type C object requires that 
much of the gas of the system has been vacated, and 
this is likely to prevent pop III star formation from pro- 
gressing efhciently. Hence, it seems unlikely that type A 
galaxies can retain their unique UV-slope signatures for 
much longer than ~ 10 Myr. 

For reference, we have in Table. [5] included the UV 
slope /3 [fx cx A^) predicted for our different pop III, II 
and I models, derive d the way this is currently done for 
galax ies at 2; w 7 (e.g. lBouwens et aT1l201Cll : lDunlop et al.l 
l2011h . Here, /3 is based on the Hub ble Space Tele- 
acojoe WFC3 J125 — ^^leo colour, using the lBouwens et all 
()2010[ ) calibration. At z = 7, these filters sample the rest 
frame UV at wavelengths of w 1370-1750 A and « 1760- 
2110 A, respectively. We caution, however, that the /3 
parameter do es depend on the exact wavelength range 
sampled (e.g. |.Raiter_et al. 2010b), especially for models 
involving nebular emission. Additional UV slope predic- 
tions for some of the pop HI SEDs used in this work are 
available in elect ronic format, along w ith numerous other 
quantities, from lRaiter et al.l ()2010bf ). 




1500 2000 2500 3Q00 3500 4000 4500 
^(A) 

Fig. 10. — The rest frame UV spectra of 1 Myr old, 
instantaneous- burst pop I (black line), pop II (red), pop III with a 
IKroupal i20mi ) IMF (green), pop III.2 (blue) and pop III.l (cyan) 
galaxies of type C. In all cases, we assume that the SED is entirely 
due to direct star light, i.e. that the leakage of Lyman continuum 
photons is so overwhelming that the contribution from nebular 
emission is negligible. The steepness of the spectra is seen to in- 
crease when going from high to low metallicities (from black to red 
and green lines), and when going from a Kroupa (200l) IMF to 
top-heavy ones (from green to blue and cyan lines). This suggests 
that pop III galaxies should be identifiable using broadband filters 
that measure the restframe UV slope. All spectra have been nor- 
malized to the same fiux at around 1300 A to allow differences in 
the spectral slope to be more easily spotted. 

III.l systems only last for about « 3 Myr before fading 
away and Pop HI. 2 systems will eventually (at ages > 10*" 
yr) attain colours redder than the W2nn — TO444 < —1.2 
criterion. Pop HI galaxies with a IKroupal (j200ll ) IMF 
evolve redward of this limit after about 4 Myr. More 
extended star formation scenarios would in principle al- 
low pop HI galaxies to retain their m2oo — TO444 < — 1.2 
colours for longer periods of time. In the case of a con- 
stant s tar formation rate, type C Pop HI galaxies with a 
iKroupa (,2001,1 IMF would display m2oo — TO444 < — 1.2 



6.1. The redshift evolution of type C galaxy colours 

Since the type C SEDs displayed in Fig. [TOl show some 
deviations from pure power-law spectra, colours such as 
'71200 ~ TO444 will evolve as a function of redshift. In 
Fig. lllb . we show the redshift dependence of ?7i2oo — '^1444 
for 1 Myr old galaxies of both type A and C. Throughout 
the redshift range z « 7-13 (shaded region), the region 
771200 ~ 77*444 £ —1.2 should effectively separate pop HI 
objects from all other types of galaxies. All type A galax- 
ies (dashed lines) appear at 771200 ~ 771444 > —0.6 within 
this redshift range. The redshift limits over which the 
777200 ~ 771444 ^ — 1-2 criterion applies are set by the fact 
that, at z < 7, pop I and pop II of type C start to oc- 
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Fig. 11. — The prospects of using the JWST/NIRCam m2oo ~ m444 colour to distinguish type C, pop III galaxies from other objects, 
a) The m2oo ~ m444 evolution of instantaneous-burst, type C populations as a function of age at 2 = 8 . The line colours indicate different 
metallicity and IMF combinations: Pop III.l (cyan line), Pop III. 2 (blue). Pop III with a IKroupal 1)20011 ) IMF (green), Pop II (red). 
Pop I (black). A Kroupa (2001) IMF has been assumed for all of the Z = 0.0004-0.020 populations. The shaded region indicates the 
'T1200 — r7i444 = — 1.2 limit below which only pop III galaxies (green, blue and cyan lines) are likely to be found. All type A galaxies appear 
at m2oo ~ m444 > —0.4 at this redshift. The error bars on the observed colours would be 0.15 mag in the case of lOo" detections in the 
relevant filters, b) The redshift evolution of the m2oo — m444 colours for 1 Myr old populations. The solid lines are the same as in the left 
panel, whereas the dashed lines represent a selection of type A galaxies; pop III. 2 (blue dashed) and pop I (black dashed). Throughout the 
redshift range z 7-13 (shaded region), the criterion m2oo —1^444 < —1.2 (thin horizontal dashed line) should effectively separate pop III 
objects from all other types of galaxies. 



cupy this region of colour space, and at z > 13, Lya 
absorption in the intergalactic medium is expected to 
cause substantial flux losses in the F200W filter, thereby 
rendering the TO200 ~ TO444 colour extremely red for all 
populations. Dust extinction, if effective in pop III galax- 
ies, would of course make the m2oo ~ mm colour red- 
der and could shroud some bona fide pop III systems 
in the colours of more metal-enriched systems. While 
it is possible that some pop III galaxies could masquer- 
ade as pop II/I galaxies (either due to old age or due to 
dust), the opposite is not true - if objects are detected 
at 771200 — TO444 < —1.2, they would be very good pop III 
candidates. 

6.2. Lingering nebular emission 

While measurements of the UV slope of type C galaxies 
could in principle help constraint the pop III IMF, the 
question remains whether pure type C, pop III galax- 
ies are at all likely to be detectable. As explained in 
Sect. [HI it would take a tremendous amount of leakage 
to bring nebular emission down to negligible levels in the 
SED of a pop III galaxy. Due to the very red colours 
of young type A galaxies (m2oo — 771444 > —0.4) and the 
high ratio of nebular to stellar flux in the F444W filter 
(Fig. U), even small amounts of lingering nebular light 
would jeopardize the possibility of identifying pop III, 
type C objects. We find that, depending on the exact 
redshift {z « 7-13), IMF and age, a Lyman continuum 
escape fraction of /esc ^ 0.95-0.99 would be required 
for a pop III galaxy to meet the m2oo — 771444 < —1.2 
limit. With such small nebular contributions to the over- 
all flux, a type C pop II galaxy would need to be an order 
of magnitude more massive than a type A object to be 
detectable with JWST (Fig. SD. While the relative con- 
tribution from nebular emission decreases at high ages 
for instantaneous burst models (> 10 Myr for IMFs that 
allow populations to remain luminous for that long) , the 



intrinsic UV slope at the same time becomes redder (see 
Fig. [TTk). thereby jeopardizing any chances of detect- 
ing a unique pop III signature. Adopting a more ex- 
tended star formation history would make pop III galax- 
ies retain their very blue UV slopes for longer, but also 
keeps the /nob//stars ratio up, and therefore still requires 
a very high fcsc to allow the broadband criteria to apply. 
Switching to colours based on JWST filters at shorter 
central wavelengths would reduce the relative contribu- 
tion from nebular emission somewhat (this is evident 
from Fig. [T]), but also decreases the difference in mag- 
nitudes between the intrinsic colours of pop III and pop 
II/I. For instance, pop III galaxies at 2; = 8 could dis- 
play 7)7150 — 771200 colours (based on fluxes in the NIRCam 
F150W and F200W filters at 1.5 and 2.0 /xm, repspec- 
tively) redder than those of pop I and II galaxies for 
~ 10*" yr even in the case of /esc ~ 0.9 (i.e. lower than 
the /esc ^ 0.95-0.99 quoted above, albeit not by much), 
but the colour difference between these models is then 
< 0.05 mag, which would be very difficult to measure in 
practice. 

Whether pop III objects with Lyman-continuum es- 
cape fractions /esc ~ 0.9-0. 99 exist is an open ques- 
tion. Recent simulations by iJohnson et al.l ([2009) sug- 
gest that /esc in this range may be produced if the pop 
III IMF is extremely top-heavy, but taken at face value, 
the prospects of identifying pop III galaxies through the 
colour signatures of type A objects appear more promis- 
ing, since they are likely to hold for a wider range of 
escape fractions (/esc ~ 0-0.5). 

7. DISCUSSION 

7.1. Caveats related to the [OIII] method 

In Sect.ini we argued that pop III galaxies with SEDs 
dominated by nebular emission (type A) could poten- 
tially be identified based on their JWST colours, since a 
combination of strong hydrogen emission lines like Ha, 
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yet absent metal emission lines like [OIIIJA5007 would 
gives rise to very peculiar colours over certain redshift 
intervals. Here, we discuss a couple of potential caveats 
with this method. 

7.1.1. The oxygen abundances of pop III galaxies 

Our tests indicate that pop III, type A signatures based 
on the lack of [OIIIJA5007 emission are erased as soon as 
the oxygen-to-hydrogen abundance relative to that of the 
Sun reaches [0/H] > -2.2 (or (0/H) > 5 x lO'^ in ab- 
solute numbers). Depending on the metallicity threshold 
at which gas switches from the formation of pop III stars 
with a top-heavy IMF to pop II with an IMF more typ- 
ical of that in the local Universe (often assumed to hap- 
pen at Z = 10"'' to 10-6, e.g. [Schneider et al.|[20Q6bl ). 
and on the oxygen yields of the first supernovae, this 
oxygen criterion may translate into a metallicity crite- 
rion for type A, pop III galaxy signatures that is either 
slightly lower or higher than that which governs the pop 
III-II t ransiti on. Hence, the method proposed by us and 
llnoud (|2011bl ) for detecting pop HI, type A galaxies may 
in principle either miss a substantial fraction of galaxies 
that are still able to form pop HI stars, or to catch galax- 
ies that have already switched to pop II star formation. 
Clearly, a method of this type is only meant to generate 
pop HI galaxy candidates. Follow-up spectroscopy will 
be required to further probe the exact nature of these 
targets. For objects that are already close to the detec- 
tion threshold for imaging in ultra deep fields, this may 
of course be very difficult to achieve with JWST itself. 
However, such follow-up observations may for instance 
be feasible with the upcoming 42 m European Extremely 
Large TelescopJ^. 

7.1.2. Ionization levels 

Another caveat associated with weak or absent 
[OIIIJA5007 emission as a criterion for pop HI galaxies 
is that pop II/I galaxies with very low ionization pa- 
rameters may also display very low [OIII] fluxes. In the 
framework of a single spherical nebulae surrounding the 
star-forming core of the galaxy, this could correspond 
to the situation where stellar feedback has increased the 
inner radius of the cloud to large radii. However, such 
objects may no longer qualify as type A objects, but 
may rather belong to type B (see Fig. [3|), with a drasti- 
cally lower nebular fluxes and potentially distinct spec- 
tral characteristics as a result. A detailed investigation 
of the properties of type B galaxies would be required to 
settle this issue. 

7.2. Narrowband imaging 

One may, perhaps, argue that narrowband surveys are 
more likely to be successful in isolating pop HI galaxies 
than any broadband survey. Aside from "classical" nar- 
rowband features of pop HI objects, like very high Lya 
or [HeII]A1640 equivalent widths, one may for instance 
look f or the Balme r jump at 0.36 /im (e.g. Schaerer: 2002, 
120031 : IInouell201 1bl. All objects with spectra dominated 
by nebular emission should display a sudden jump in the 
continuum level at the short-wavelength side compared 
to the long wavelength side of this limit. Since metal- 
enriched stars are less efficient in ionizing hydrogen and 

http:/ /www. eso.org/sci/facilities/eelt/ 



have redder stellar continua, they will - relatively speak- 
ing - be less dominated by nebular emission (see Fig. [2]) . 
Therefore, one may conjecture that the objects with the 
most pronounced Balmer discontinuities are likely to be 
pop HI galaxies, and that suitably placed narrowband 
filters would be able to pick up on this. However, our 
models indicate that the difference between the Balmer 
jumps between pop II and pop HI galaxies are so small 
that it will be next to impossible to tell them apart using 
this feature. While the Balmer jump may still be used for 
identifying extremely young, nebular-dominated objects 
at low metallicities (albeit not necessarily just pop HI 
objects), the redshift range over which this feature can 
be efficiently probed with a given setup of JWST nar- 
rowband filters is very small (Az ~ 0.1). This makes the 
number of potential targets in the corresponding volume 
so tiny, that it would seem more time-efficient to use 
broadband criteria to pre-select pop HI galaxy candi- 
dates and then to go directly for follow-up spectroscopy. 

7.3. How much do the details of the pop III IMF 
matter? 

iBromm et ahl (|200l argue that very massive pop HI 
stars are self-similar, in the sense that a fixed mass of 
gas converted into either 100 M©, 300 M© or 500 Mq 
stars would produce nearly identical SEDs. Indeed, the 
lifetimes, ZAMS temperatures and the time-averaged, 
mass-normalized Lyman continuum fluxes of M sa 100- 
5OO_M0 stars differ by no more than 10-30 % l|Schaereii 
l2002f ). However, the approximation of self-similarity de- 
teriorates once pop HI stars at M < 100 Mq are consid- 
ered, and does not hold for the Hell ionizing spectrum 
and the corresponding recombination lines, or for stars 
that have evolved off the main sequence, as discussed in 
detafl by .Schaerer, (.200^ . 

As seen in Fig. |31 the mass-to-light ratios produced 
by our pop III.l (~ 100 Mq) and pop III.2 (~ 10 Mq) 
IMFs differ by a factor of « 2-3. This is mostly due to 
the inclusion of low-mass stars which contribute substan- 
tially to the mass and not to the luminosity. Hence, the 
colours produced by the different pop HI IMFs remain 
very similar fFigs.l6l[7l[8land [TT|) . For a flxed star forma- 
tion history, galaxies with pop HI. 2 IMFs can continue 
to shine for longer than those with pop III.l IMFs, but 
their colours are ot herwise identical. Pop HI galaxies 
with ' Kroupal (j2001| ) IMFs display slightly less extreme 
colours, but this is largely caused by a shift in the time 
scales - pop HI. 2 galaxies will display the same colours, 
just at a slightly higher age. It is only in the case of a 
purely stellar SED (our case C; Fig [TT]) that Pop III.l 
and pop HI. 2 get separated from the lKroupal (j2001[ ) IMF 
even at very low ages. Hence, the detection of objects 
with extreme colours would support the idea of a pop HI 
IMF wi t h mor e high-mass stars than predicted by the 
iKroupal (|200lD IMF, but the prospects of disentangling 
a pop III.l IMF from a pop HI. 2 IMF using integrated 
colours alone seem bleak. 

8. SUMMARY 

Yggdrasil is a new spectral synthesis model which in- 
cludes pop I, II and HI stars, dark stars, nebular emis- 
sion and dust extinction. Using this model, we derive 
the stellar population masses of the faintest galaxies de- 
tectable by JWST through broadband imaging in ultra 
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deep fields. Assuming 100 h exposures per filter, we find 
that JWST may be able to detect pop III galaxies with 
stellar masses as low as ^ 10^ M© at z = 10, whereas the 
corresponding limit for pop II/I galaxies is ~ lO^M©. We 
also argue that the broadband fluxes of young (age < 10 
Myr) pop III galaxies are likely to be significantly af- 
fected by nebular emission, unless the fraction of ionizing 
radiation escaping directly into the intergalactic medium 
is extremely high (/esc ^ 0.95). 

We also discuss different strategies for selecting Pop 
III galaxy candidates based on their JWST broadband 
colours, both in the limiting cases of having dominant 
(low /esc) and negligible (very high /esc) nebular emis- 
sion. In the former case, we particularly highlight the 
possibility of adding imaging in two MIRI filters (F560W 
and F770W) to the NIRCam filter sets usually discussed 
in the context of JWST ultra deep fields, since this should 
allow for a clean selection of pop III galaxies at z ss 7-8. 
Selecting targets for spectroscopic follow-up using colour 



criteria based solely on NIRCam filters would be more 
economical in terms of exposure time, but also increases 
the risk of interlopers. In the case of pop III galaxies 
dominated by direct star light (very high /esc), we argue 
that imaging in the NIRCam F200W and F444W filters 
should allow pop III candidates to be selected through- 
out the redshift range z w 6.5-13. 
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